The Mott insulator transition in a two dimensional atomic Bose gas 
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Cold atoms in periodic potentials are remarkably versatile quantum systems for implementing 
simple models prevalent in condensed matter theory. Here we realize the 2D Bose-Hubbard model 
by loading a Bose-Einstein condensate into an optical lattice, and we study the resulting Mott 
insulating state (a phase of matter in which atoms are localized on specific lattice sites) . We measure 
momentum distributions which agree quantitatively with theory (no adjustable parameters). In 
these systems, the Mott insulator forms in a spatially discrete shell structure which we probe by 
focusing on correlations in atom shot noise. These correlations show a marked dependence on 
the lattice depth, consistent with the changing size of the insulating shell expected from simple 
arguments. 
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In recent years, ultra-cold atoms confined in opti- 
cal lattices have realized fascinating strongly interacting 
condensed matterphenomena, including the Girardeau- 
Tonks gas in ID [Ji and the superfluid to Mott insula- 
tor transition in 3D, and ID 0, The Mott-insulating 
state is a strongly correlated phase of matter where in- 
teractions localize the constituent particles to individual 
lattice sites. While the existence of a 2D Mott insu- 
lator has been verified in the cold atom system it 
has gone largely unexplored. In all of these systems, an 
atomic Bose-Einstein condensate (BEG) is loaded into 
the optical lattice formed by interfering laser beams, con- 
fining the atoms in a nearly perfect 1, 2, or 3D peri- 
odic potential. The atom-optical system differs substan- 
tially from traditional condensed matter systems both in 
control (the trapping potentials can be changed dynam- 
ically on all relevant experimental time scales) and in 
measurement opportunities (e.g., imaging after time-of- 
flight provides a direct measurement of the momentum 
distribution). Imaging complements techniques usually 
available in condensed matter systems, and allows us to 
study the 2D Mott insulator via its momentum distribu- 
tion and correlations in its noise @, 0, [1, • Our mea- 
sured momentum distributions agree quantitatively with 
the predictions of a 2D theory with no adjustable param- 
eters (even approaching the insulator-superfluid transi- 
tion, where perturbation theory breaks down). In the 
cold atom system, the transition from superfluid to insu- 
lator occurs smoothly; the system segregates into shell- 
like domain(s) of insulator and superfluid 13, U, 12, T3 |. 
Even in the absence of direct imaging, we flnd signatures 
of the insulating shell in noise-correlations. These noise- 
correlations dependence on the lattice depth, trending as 
expected due to the changing size of the Mott insulator 
region. 

Our sample consists of an ensemble of 2D Bose systems 
in a combined sinusoidal plus harmonic potential. Absent 
the harmonic confinement, the system is well described 
by the homogeneous Bose-Hubbard Hamiltonian l3, 3] , 
which models bosons on a lattice in terms oft, the matrix 



element for tunneling between sites, and U , the on-site 
energy cost for double occupancy. At sufficiently low 
temperature, this system is either a superfluid (SF) or a 
Mott insulator (MI). In 2D, the SF state has quasi- long 
range phase coherence and a gapless excitation spectrum 
linear in momentum; in contrast, the MI phase has a gap 
between the ground and flrst excited states. For unit fill- 
ing (on average one atom per lattice site), the 2D Bose- 
Hubbard system exhibits a zero temperature quantum 
phase transition from SF to MI |U] when the dimen- 
sionless ratio t/U « 0.06 [H, which we denote by 
{t/U)c- When t/U vanishes, the unit filled MI has ex- 
actly one atom per-lattice site. At finite t/U this is not 
fixed and the variations from unit occupation are corre- 
lated between nearby sites. This is associated with coher- 
ence between these sites, resulting in broad, diffractive, 
structure in the momentum distribution 



In current experiments, the additional harmonic trap- 
ping potential modifies the physics from the ideal case 
described by the homogeneous Bose-Hubbard Hamilto- 
nian. In the trapped system, domains of the insulating 
phase grow continuously as t/U decreases (we vary t/U 
by changing the depth of the lattice potential), forming 
a spatially discrete shell structure [l3, [TH, [H, [3] ■ This 
is most simply interpreted in terms of a "local chemical 
potential" which varies from a peak value in the cen- 
ter of the trap to zero at the edges (i.e., a local den- 
sity approximation, LDA). In this case, the system seg- 
regates into one or more shells of MI separated by SF. In 
our experiment, the coexisting phases are SF and unit- 
occupancy MI. Examples relevant to our experiment are 
shown schematically in the top row of Fig. [1] where (A) 
is a system which has a SF core enveloped by a MI shell. 
(B) and (G) represent deeper lattices where the sample 
is MI nearly throughout. 

We prepare the Bose-Hubbard system at a specific 
value of t/U (28|, measure the momentum distribution, 
and extract correlations in its noise. The characteristic 
loss of diffraction in the momentum distribution as 
the system goes deeper into the Mott regime is shown in 
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FIG. 1: Results for 3 different values of t/U. Top row: expected in situ density profile from a LDA calculation using the 2D MI 
phase diagram from ref. [1^. The regions of dark grey denote the MI phase, and the light regions indicate SF. In all cases, the 
MI is surrounded by a small ring of SF. Middle row: imaged atom density versus momentum. Bottom row: noise correlations 
of the images as a function of momentum difference. Each displayed image was averaged from about 60 raw images; also, to 
reduce the noise, the noise-correlation data shown was averaged with itself, rotated by 90° . While the momentum distributions 
at t/U = 4 X 10""^ and 0.5 x 10"'' are quite different, the shell structures hardly differ. 



middle row of Fig. [TJ At low temperature, the MI is ex- 
pected to form in a shell structure, with each region con- 
tributing to the image; it can be difficult to deconvolve 
the separate contributions, particularly near the transi- 
tion [Si- Despite these limitations, we find that a 
theory for the homogeneous system is surprisingly good 
at describing the full momentum distribution. 

The bottom row of Fig. [1] shows noise correlations 
of the momentum distribution, whose diffractive nature 



does not vanish, even deep in the Mott regime [gI, ^ . The 
width and area of the diffractive noise correlation peaks 
should depend on the size of the Mott region. Deep 
in the Mott limit, the width of the correlation peaks 
should be diffraction limited by the inverse linear size of 
the Mott domain, and the area of the correlation peaks 
should depend on the number of atoms in the MI. Indeed, 
we observe that the area of the correlation peaks clearly 
increases with increasing t/U, while the width shows a 
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more subtle trend. 



EXPERIMENTAL PROCEDURE 



^■^Rb BECs with Nt_ 
l,mp = —1) state \1{ 



We produce nearly pure 3D 
1.7(5) X 10^ ^ atoms in the : 
The 3D BEC is separated in 200 ms into an array of 
about 60 2D systems by an optical lattice aligned along 
z (vertical lattice) . This lattice is formed by a pair of lin- 
early polarized A = 820 nm laser beams [3(|. In addition, 
a square 2D lattice in the x-y plane is produced by a sec- 
ond beam arranged in a folded-retroreflected configura- 
tion [2o| ; the polarization is aligned in the x-y plane. This 
x-y lattice is gradually applied in 100 ms 3l|. The inten- 
sities of the vertical and x-y lattices follow exponentially 
increasing ramps (with 50 ms and 25 ms time constants 
respectively) which reach their peak values concurrently. 
These time-scales are chosen to be adiabatic with respect 
to mean-field interactions, vibrational excitations, and 
tunneling within each 2D system . Ideally, the loading 
process transfers a ground state BEC to the ground state 
of the lattice system. (While we do not know the temper- 
ature in the lattice, we begin with a nearly pure BEC and 
are confident that our loading procedure does not cause 
excessive heating. To check this adiabaticity, we created 
a MI as described below, then decreased the lattice po- 
tential in about 30 ms and verified the reappearance of 
the sharp diffraction orders indicative of the SF phase.) 
The final vertical lattice depth is always 30(2)i5/j; the 
single photon recoil momentum and energy are defined 
as kR = 27r/A and Er = h^kn^ /2m = /i x 3.4 kHz. The 
final depth of the x-y lattice determines t/U and ranges 
from V = Q to 31(2)i?fl. In our experiments, the lattice 
depths are measured by pulsing the lattice for 3 ^s and 
observing the resulting atom diffraction 21 1 . 

Once both lattices are at their final intensity, the 
atomic system consists of a set of 2D gasses each in a 
square lattice of depth V and with a typical density of 
1 atom per lattice site. The atoms are held for 30 ms 
at which point all confining potentials are abruptly re- 
moved (the optical lattice and magnetic potentials turn 
off in < 1 /iS and ~ 300 /iS, respectively). When the 
potentials are thus "snapped off" the initially confined 
states are projected onto free particle states which ex- 
pand for a 20 to 30 ms time-of-flight (TOF); they are 
then detected by resonant absorption imaging (see the 
supporting Methods section for details). Because initial 
momentum maps into final position, each image corre- 
sponds to a measurement of the 2D momentum distribu- 
tion n{kx, ky). Ideally, when averaged over many realiza- 
tions, this converges to the expectation value {fik^.ky} of 
the 2D density operator. 

The middle row of Fig. [T] shows a series of such av- 
erages, at different lattice depths, starting near the MI 
transition and crossing deep into the MI phase 3] . These 



images illustrate the smooth progression from sharp 
diffraction peaks on a small background (A), to broad 
diffraction peaks with a considerable background (B), 
culminating with a near-perfect Gaussian distribution 
(C). Because the system is divided into domains of SF 
and MI, the diffraction images can be difficult to inter- 
pret - only returning to simplicity in the limits of a shal- 
low or deep lattice when the entire system consists of 
just one phase (SF or MI). Understanding quantitatively 
the intermediate regime requires careful modeling of the 
inhomogeneous density distribution [l8| . 



MOMENTUM DISTRIBUTIONS 

In the deep lattice limit, the MI wave function |\l/o) 
has an exact number of atoms in each occupied site 
(a Fock state) - in our case n — 1. At large but fi- 
nite depth, the unit-occupied state is modified to first 
order in perturbation theory with a small mixture of 
neighboring particle-hole pairs. As in a 3D lattice, this 
gives a modulated momentum distribution, here (hk) = 
N\w{k)f{l-^a[coa{7rk^/kr)-^coa{nky/kr)]} [13, \M- 
Here a — 8t/U, and it;(k) is the Fourier transform of the 
localized Wannier states in the optical lattice sites; ^(k), 
determined by the lattice depth, is well approximated by 
a Gaussian for all data described here. Because the den- 
sity modulations result from an interference between the 
unit-occupied Mott state and the particle-hole admix- 
ture, the lowest order correction to (uk) is proportional 
to t/U, even though the probability for double occupancy 
scales as {t/U)"^. 

This modulation in density has been verified to be first 
order in t/U over a range of parameters in 3D 



17|. As 



with the 2D case, which we discuss below, this agreement 
near the MI-SF transition can be surprising for two rea- 
sons: (1) as seen in the top row of Fig. [1] only a fraction 
of the inhomogeneous system may be in the MI phase; (2) 
as t/U approaches the critical value from below, higher 
order contributions become important. 

To quantify the next order term, 
analytic result (in the random phase 
22| ) to second order in t/U. 
order term in the momentum 
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we expand an 
approximation. 
The expected 
distribution is 
This yields ad- 
- cos{2'nky / kr)) , 



72{t/U)'^[cOs{TTka:/kr) + COs{TTky / kr)f 

ditional Fourier terms Pi{cos{2'!rkx / kr 
and P2C0s{'Kkx/kr) cos^nkx/kr)] we define the average 
coefficient p ^ {Pi + /32)/2 = 90{t/U)^. (Terms higher 
order in t/U also contribute slightly to these Fourier com- 
ponents.) 

We extract the Fourier coefficients a and /? from in- 
dividual images and plot them in Fig. [2l The top panel 
shows a versus t/U; the red dashed line is the theoreti- 
cal prediction with no adjustable parameters: a = 8t/U. 
The vertical black line denotes {t/U)c = 0.06, where we 
expect the MI to first appear somewhere in our sample. 
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FIG. 2: Top: the symbols show the measured first order coef- 
ficient a versus t/U. For reference the red dashed line shows 
the expected value a = 8t/U. Each point is acquired from 
a single image at known lattice depth, therefore known t/U. 
(A), (B), and (C) indicate the value of t/U for the three sets 
of data in Fig. [T] Bottom: the symbols denote the averaged 
second order coefficient (3 versus t/U. The red dashed line 
is the expected power-law /3 = QOjt/U)^. In both figures the 
vertical dotted line at {t/U)c = 0.06 0,131, indicates the ex- 
pected location of the 2D SF-MI transition. (Some extracted 
values of /3 are negative, and to not appear on the log plot.) 



Points to the left of the black dotted hne (MI regime) 
agree with perturbation theory. 

The bottom panel of Fig. [2] shows the /3 versus t/U. 
The red dashed line is the predicted value /3 = 90{t/U)^. 
In the MI regime, our measurements agree with expec- 
tations, within experimental uncertainty. Note that near 
the MI transition the first and second order term be- 
come comparable, indicating the incipient breakdown of 
perturbation theory. 



t/U =25x10" 
V =12.5E. 




0.0 
0.3 

0.2 

0.1 

0.0 



t/U =7x10 
V =17Er 





-1.0 -0.5 0.0 0.5 
Momentum [/c^] 



1.0 



FIG. 3: Cross sections of normalized quasi-momentum distri- 
butions (along x + y) at three values of t/U (along this diago- 
nal the section extends to = ±\/2). In each case the data is 
plotted along with the theoretical profile (red lines) [23] . The 
dashed lines (not visible in the bottom panel) refiect the un- 
certainty in the theory resulting from the single-shot ±0.5Er 
uncertainty in the lattice depth. 



In Fig. [3l we directly compare the measured quasi- 
momentum distribution with theory for three values of 
t/U . The modulated momentum distribution is the prod- 
uct of the magnitude squared Wannier function |w(/c)p 
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and the quasi-momentum distribution (periodic along x 
and y in the reciprocal lattice vectors, 2kfi). To ex- 
tract the normalized quasi-momentum distribution we 
divide the data by [^(/c)^, then calculate a properly 
weighted average of all points separated by multiples of 
the reciprocal lattice vectors, and finally normalize. The 
data in Fig. [3] are cross sections of the quasi-momentum 
distribution along x + y plotted with the prediction of 
the RPA theory (red lines). At low t/U the quasi- 
momentum distributions are cosinusoidal. As t/U in- 
creases toward (t/U)c, contributions of higher Fourier 
terms become important, as is evident in the cross sec- 
tion at t/U = 25 X 10~^ which deviates from a cosine. 
The shape of the measured distribution matches the pre- 
dictions of theory with no free parameters. 



NOISE CORRELATIONS 

Although our sample is inhomogeneous, a 
homogeneous-system theory provided a remarkably 
good representation of the data discussed above. This 
results from two facts: (1) except quite close to the 
MI transition, nearly all of the system should be unit 
occupied MI, and (2) by focusing only on the largest 
momentum scales (corresponding to spatial length 
scales on the order of one- or two- lattice sites) we are 
insensitive to the size of the MI. In contrast, the size of 
the MI almost exclusively determines the area and width 
of the peaks in the noise correlations signal (bottom row 

of Fig. [ID gg. 

We determine noise correlations from our images of 
atom density n{kx, ky) by computing the autocorrelation 
function averaged over many images: 



S{dkx,6ky) 



-'0 



^(^a^; ky^Tii^k'j; Sk^^^ ky — {~ Sky^dk'^dkj 



we 
the 



normalize by the autocorrelation of 
average {n{k^,ky)), so So{Skx,Sky) = 



J {n{kx, ky)) {n{kx + Skx, ky + Sky)) dkxdky. In other 
words, we determine S {Skx, Sky) by calculating the au- 
tocorrelation function (ACF) of each image separately, 
and then average over many realizations, typically (40 
to 80). (We normalize this by the ACF of the averaged 
images to remove the dependence of S{Skx,Sky) on the 
momentum distribution {n{kx,ky)).) It can be shown 
that this quantity has diffractive structure when |\E'o) is 
a pure Fock state [1, [i^], with noise correlation peaks 
separated by 2kii (the reciprocal lattice vector). The 
spacing of the peaks therefore reveals the underlying 
lattice structure, even in the absence of diffraction in 
the momentum distribution. In addition, the correlation 
peak areas A and widths Sk, provide information about 
the system. 

In the limit of a deep lattice, the ground state of 
our system is a 3D array of lattice sites with exactly 



one atom per site. As is usual for diffraction phenom- 
ena, the width of the correlation peaks is determined by 
the size of the array, and is proportional to L~^; where 
L oc N^^^ is the linear extent of the MI region and N is 
the number of sources (lattice sites). Likewise, the area 
under a peak is related to the atom number in the MI 
hy A = {2kf()'^/N Q. As with most noise, the noise- 
correlation variance (the square of the standard devia- 
tion), scales as N. We normalize by a quantity which 
scales like 1/iV^, so the relative fluctuations described 
by A have an overall 1/N dependence. We calculate 
that this remains true including the order t/U correc- 
tion to the MI state, so even to first order, noise cor- 
relations contain information only about the size of the 
Mott-insulating system. (Order (t/U)'^ corrections to the 
noise-correlations, which we have not investigated, may 
alter this behavior). 

As the system approaches the critical point from the 
MI side, we expect the size of the MI region to shrink, and 
correspondingly the correlation peaks to become wider, 
with increasing area (the total number of atoms in the 
experiment Nt remains fixed; only the fraction of atoms 
in the MI N/Nt decreases). Figure |4] shows this general 
behavior. Our study of the noise correlations near the MI 
critical point is enabled by a masking procedure: in each 
image we eliminate regions of radius 35 /Ltm « 0.3 x ka 
centered on the diffraction peaks before computing cor- 
relation functions. This removes spurious effects of the 
sharp diffraction peaks evident near the MI-SF transition 
(see supporting Methods section). 

Figure [4lA_ shows the area under the noise correlation 
peaks measured by fitting the peaks to 2D Gaussians. 
The data show that as t/U increases the area also in- 
creases; this indicates that the fraction of the system in 
\the MI state is decreasing as expected. The area ex- 
, PQcted for small t/U is {2kii)^ /Nt, however, our data 
xends to about 45% of this in a deep lattice (a similar 
suppression of the noise signal was observed in Ref. f^). 
We attribute at least some of this discrepancy to colli- 
sions during the ballistic expansion of the system, which 
modify some atoms' trajectories, removing those atoms 
from the correlation features; errors in number calibra- 
tion could also contribute. 

The data are plotted along with a red dashed line show- 
ing the change in peak area expected if it were due solely 
to the finite size of the Mott domain, which we calculated 
in a LDA using the 2D MI phase diagram from ref. [l5| . 
After we scale the model by a factor of 0.45 to account for 
the overall decrease in the measured correlation-signal, it 
agrees to within our uncertainties. 

The width of the correlation peaks are shown in Fig.|4(3 
where the symbols are the measured RMS peak-widths 
from a Gaussian fit, and the red dashed line is the ex- 
pected peak-width for a pure MI with size given by our 
LDA model. At small t/U, the experimental data sat- 
urate to about 0.045/cfl, compared with the 0.015/cfl ex- 
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FIG. 4: Measured average area and width of the two central 
independent noise-correlation peaks. The uncertainty bars re- 
flect the statistical uncertainty of the fit due to background 
noise. Top: average area expressed in units of kn^ /Nt- In 
these units, the expected area for a perfect n — 1 Fock state is 
4; the data tend to about 1.8 kn^ /Nt in a deep lattice. The 
red dashed line denotes the expected area due to the finite 
number of atoms in the MI as calculated in our LDA model, 
which was scaled by 0.45 to lie upon the data. Bottom: the 
solid symbols are the measured width in ka of the correlation 
peaks. The red dashed line is the expected peak-width based 
on our LDA computation of the system size. The black dashed 
line shows the expected peak- width including the imaging res- 
olution of O.OSfefl (see supporting Methods section). In both 
panels, the vertical dotted line show {t/U)c ~ 0.06 [isl. [l^. 
the expected location of the 2D SF-MI transition. 



these effects would separately limit the measured peak- 
width to about 0.03 and 0.04 fc_R, respectively. (The 
width may also be influenced by mean field during ex- 
pansion.) The black dashed line shows the modeled peak- 
width added in quadrature with a width of 0.05 kji, the 
quadrature sum of the effects described above. The width 
of the correlation peak appears to increase more rapidly 
than expected from size arguments alone. This suggests 
that the peak-width may depend on the detailed proper- 
ties of the MI state (as the momentum distributions do). 
It is also possible the the finite temperature of our sys- 
tem may influence the peak- widths, or that our masking 
procedure may introduce an unknown systematic error 
in the peak-width (in spite of evidence to the contrary: 
see the supporting Methods section). 

In our simple calculations (valid to first order in t/U), 
the area and width of the correlation peaks depend only 
on the size of the Mott domain, computed in a LDA. 
More sophisticated theoretical techniques can be applied 
to this problem, indeed explicit numerical calculations for 
a harmonically confined 2D Bose-Hubbard model verify 
that the peak- width increases as t/U approaches the MI- 
SF critical point from the MI side 2J]. Additionally, 
these authors find that when the SF phase first appears 
in the core of the MI (as in our experiment) the peak- 
width increases more rapidly as t/U approaches the crit- 
ical value [l^]. Numerical results for ID systems also 
show a increase in the area dst/U increases to the critical 
point [2^, [2^ . The parallels between our measurements 
and these numerical results indicate the need for further 
study of noise correlations near the Ml-transition. 

In this paper we have demonstrated a remarkable 
agreement between experiment and theory describing the 
momentum distribution in a 2D MI over a wide range of 
conditions, and to second order in perturbation theory. 
Additionally, we see that correlations in the atom shot 
noise can be a tool for probing the SF-MI phase tran- 
sition, and yield information about the fraction of the 
system in the MI. Even when the momentum distribu- 
tion is featureless, the noise-correlations show the lattice 
structure and indicate system size. This adds support 
to proposals to identify the phases of extended Bose- 
Hubbard models (including a possible supersolid phase), 
using a combination of momentum and noise-correlation 
measurements ^24]. 

We would hke to thank V. W. Scarola, E. Demler, A. 
M. Rey, and C. Williams for conversations, M. Anderlini 
and J. Sebby-Strabley for lattice setup. We acknowledge 
the financial support of ARDA, and ONR; and I.B.S. 
thanks the NIST/NRC fellowship program. 



pected in our model. This saturation is due to at least 
two effects: (1) the finite resolution of our optical system 
(see Methods supporting section), and (2) the ~ 15 /xm 
initial radius of the sample. We estimate that each of 
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METHODS 
Image acquisition 

Once the confining potentials are removed the atom 
cloud expands for a TOF (20.1 ms for the data prob- 
ing for noise-correlations, and 29.1 ms otherwise). The 
atoms are then illuminated for tpuisc = 50 /zs by 
a circularly polarized probe beam resonant with the 
\F — 2,mF = 2) \F = 3,mF = 3) cycling transition of 
the ^^Rb D2 line. The probe beam was carefully aligned 
parallel with the vertical lattice, and approximately with 
a « 0.2 mT bias field, with intensity / = 0.2/o, where Iq 
is the resonant saturation intensity. During the imag- 
ing pulse an atom with lifetime F will scatter iV^ = 
rtpuisc/[2(l + /o/-^)] ~ 160 photons out of the beam. We 
image the shadow cast by the atoms onto a CCD array. 

Our ultimate imaging resolution is set by several ef- 
fects: the resolution limit of the optical system (diffrac- 
tion and aberration), the finite size of the CCD pixels, 
and the movement of the atoms during imaging (which 
is negligible in our case). 

We measure a RMS optical resolution of 5 /Km by fit- 
ting a Gaussian to the spatial power spectra of atom 
shot-noise. The diffraction limit of our optical system 
is predominately determined by a pair of lenses with 
//# = 2.5. At 780 nm this pair gives a diffraction lim- 
ited spot with a 2.5 /im RMS radius. After the 6x magni- 
fication of our imaging system, the CCD has an effective 
pixel size of 2.5 /xm at the cloud. Based on the quadra- 
ture sum of these terms, a single atom should be imaged 
within a circle with a RMS radius of about 4 /xm. 

While the discussion in the body of the text focused 
upon fundamental atom shot noise, any real experiment 
is also beset with technical noise. In our absorption im- 
ages, the CCD is illuminated by a probe beam that is only 
slightly attenuated by the atoms. The CCD array counts 
the remaining photons, which arrive stochastically - pho- 
ton shot noise - which acts as technical noise. Second, 
a readout error is introduced when the photo-electrons 
in each CCD bin are counted. Together these give an 
expected RMS background statistical uncertainty in our 
reported optical depths of 0.02. In our data, a resolution- 
limited spot might contain « 25 atoms, which has a 5 
atom counting uncertainty giving rise to an optical depth 
change of 0.02. Only a fraction of that noise is correlated, 
requiring us to average many images to extract the cor- 
related noise signal. 



Image processing 

The atom noise signal represents the measured shot to 
shot fluctuations of the momentum distribution around 
the average distribution. 



The noise-correlations in the momentum distribution 
are a subtle feature which can easily be mimicked by 
correlations in the average momentum distribution if the 
appropriate average is not subtracted. Therefore, it is 
critical to accurately determine an average distribution, 
{n{kx, ky)), to subtract from the measured distributions. 
Ideally, {n{kx,ky)) is determined by averaging a large 
number of measured distributions under identical condi- 
tions. In practice, fluctuations in initial conditions (lat- 
tice depth: ±10%, atom number: ±20%, transverse ve- 
locity: ±250 /zm/s, and etc.) make it difficult to pre- 
cisely define {n{kx,ky)). As a result, spurious "ghosts" 
of diffractive structure are replicated in the correlation- 
data due to imperfect removal of the atom background. 
This problem is not evident except when sharp diffrac- 
tive peaks appear in the momentum distribution, i.e., 
t/U > {t/U),. 

To mitigate this, each image is processed as follows: we 
first fit to an empirical multi-parameter function which 
includes diffraction peaks. When subtracted from the ini- 
tial image, this removes the average atom density (ideally 
leaving only the noise signal including both atom shot- 
noise and technical noise). However, because our fitting 
function is imperfect, some fine structure still persists at 
the locations of the diffraction peaks after subtraction. 
If ignored, these features give rise to a spurious signal in 
the ACF exactly at the wave- vectors of interest. We note 
that these features exist only at the location of diffraction 
peaks and are separated by multiples oi2kR. 

We remove the effect of this localized structure on the 
ACF by masking circular regions with a 35 /im radius 
radius centered on the location of the diffraction peaks 
(the diffraction peaks have a typical radius ranging from 
15 to 20 /im). We find both empirically (when there is no 
diffractive structure, i.e., at very small t/U), and in sim- 
ulation that such masks have no effect on the normalized 
ACF. (Foiling et al used a similar technique, but only 
deep in the MI phase, to show that the correlation-peaks 
do not result from spurious effects of remnant diffrac- 
tion We verified that the size of the mask had little 
effect on the final ACF when the radius ranged from to 
45 when t/U is small (the least diffractive data) and 
from 30 to 45 /im when t/U is large (the most diffractive 
data). For masks larger than about 45 /im, we could no 
longer distinguish the desired correlation features from 
noise in the most diffractive data. 

When left unmasked, the spurious diffractive signal in- 
troduces excess correlations which (1) increase the overall 
magnitude of the correlation peaks, and (2) increase the 
width of the correlation peaks. 

There are two additional imaging artifacts which we 
remove. First we remove a high spatial frequency com- 
ponent due to noise in the readout electronics, which is 
present in every image (with variable phase and ampli- 
tude). Second, we remove a low frequency modulation 
(w 3 oscillations over the whole image) due to diffraction 
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fringes of the probe beam. 



